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Abstract 

The actin cytoskeleton in living cells has many types of crosslinkers. The mechanical interplay 
between these different crosslinker types is an open issue in cytoskeletal mechanics. We develop a 
framework to study the cooperativity and redundancy in the mechanics of filamentous networks 
with two types of crosslinkers: crosslinkers that allow free rotations of filaments and crosslinkers 
that do not. The framework consists of numerical simulations and an effective medium theory on 
a percolating triangular lattice. We find that the introduction of angle-constraining crosslinkers 
significantly lowers the filament concentrations required for these networks to attain mechanical 
integrity. This cooperative effect also enhances the stiffness of the network and suppresses non-affine 
deformations at a fixed filament concentration. We further find that semiflexible networks with only 
freely-rotating crosslinks are mechanically very similar to compositely crosslinked flexible networks 
with both networks exhibiting the same scaling behavior. We show that the network mechanics can 
either be redundant or cooperative depending on the relative energy scale of filament bending to the 
energy stored in the angle-constraining crosslinkers, and the relative concentration of crosslinkers. 
Our results may have implications for understanding the role of multiple crosslinkers even in a 
system without bundle formation or other structural motifs. 
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I. INTRODUCTION 



The mechanical response of most cells arises from the mechanics of its cytoskeleton, a 
polymeric scaffold that spans the interior of these cells, and its interaction with the extra- 
cellular environment. The cytoskeleton is made up of complex assemblies of protein filaments 
crosslinked and bundled together by a variety of accessory proteins. For example, there are 
approximately 23 distinct classes of accessory proteins such as fascin, a-actinin, and filamin 
A [l| that crosslink filamentous-actin (F-actin), a major component of the cytoskeleton 
that is resposible for the mechanical integrity and motility of cells. Given the multitude of 
crosslinkers, several natural questions arise: Are the different types of crosslinkers redundant, 
or do they each serve specific functions? Do they act independently or cooperatively? What 
are the consequences of their mechanics for the mechanical integrity and response of the cell? 
A mutation study of dictyostelium discoideum cells lacking a particular actin crosslinking 
can still grow, locomote, and develop, though with some defects, thereby suggesting at 
least partial redundancy in the crosslinker's mechanical function [2J . On the other hand, 
two types of crosslinkers working cooperatively may produce enhanced mechanical response. 
This cooperativity has been demonstrated in stress fibers crosslinked with the actin binding 
proteins (ABP) a-actinin and fascin, where stress fibers containing both a-actinin and fascin 

n 

were more mechanically stable than stress fibers containing only a -actinin or fascin |3|. 
In addition, it has been found that two different crosslinkers are required for actin bundle 

n 

formation in vivo [4|]. It could also be the case that different crosslinkers work independently 
of one another such that the dominant crosslinker dictates the mechanical response of the 
network [5] . Given these various possibilities, how the cell uses different crosslinking proteins 
to optimize for certain mechanical characteristics is an important open issue in cytoskeletal 
mechanics. 

Here, we address this redundancy versus cooperativity issue by studying a model net- 
work of semiflexible filaments crosslinked with two types of crosslinkers. We first study 
the mechanical properties of the model network with one type of crosslinker and then add 
the second type of crosslinker and look for mechanical similarities and differences with the 
original model network. In addition, we also address the redundancy versus cooperativity 
issue of two types of crosslinkers for networks made of flexible filaments. 

As for the two types of crosslinkers, we consider crosslinkers that allow the crossing fila- 
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ments to rotate freely (freely-rotating crosslinks) and crosslinkers that constrain the angle 
between two filaments. The ABP a-actinin is a candidate for the former type of crosslink- 
ing mechanics: optical trapping studies demonstrate that two filaments bound by a-actinin 
can rotate easily 6|. As an example of the latter, we consider filamin A (FLNa), which 
binds two actin filaments at a reasonably regular angle of ninety degrees, suggesting that 
FLNa constrains the angular degrees of freedom between two filaments [7f. Here, we do not 
take into account the possible unfolding of FLNa since the energy to unfold filamin A is 
large j3"9], nor do we take into account the kinetics of FLNa since we seek to understand 
fully the mechanics in the static regime first. There exist other possible examples of angle- 
constraining crosslinkers such as Arp2/3 that serves a dual role as an F- actin nucleator and 
a crosslinke r ll 0f | - While its role as a nucleator has been emphasized in lamellipodia for- 
mation fill . [12! , its role constraining the angle between the mother and daughter filaments 
is presumably also important for lamellipodia mechanics. Better understanding of the me- 
chanical role of Arp2/3 in lamellipodia may also help to distinguish between the dendritic 
nucleation model for lamellipodia formation and a new model where Arp2/3 only nucleates 
new filaments but does not produce branches fl3 |. 

In studying the mechanical properties of compositely crosslinked filamentous networks, 
we focus on the onset of mechanical rigidity as the filament concentration is increased above 



3J-|20|. Above 



some critical threshold. This onset is otherwise known as rigidity percolation 
this critical threshold, both experiments and theoretical studies of F-actin networks have 
observed distinct mechanical regimes. For dense, stiff networks the mechanical response 
is uniform or afline and the strain energy is stored predominantly in filament stretching 
modes. While for sparse, floppy networks one finds a non-affine response dominated by 
filament bending where the observed mechanical response of the network is inhomogeneous 
and highly sensitive to the lengthscale being probed 
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25| . It has been recently reported 



that there exists a bend- stretch coupled regime for intermediate crosslinking densities and 
filament stiffnesses 26]. 

While considerable progress has been made in understanding the mechanics of cytoskeletal 
networks that are crosslinked by one type of crosslinkers, compositely crosslinked networks 



are only beginning to be explored experimental! 



with one type of crosslinker theoretically 
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as are composite filament networks 



Here we investigate the mechanics of such networks as a function of the concentration 
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and elasticity of the crosslinkers and the filaments. 

II. MODEL AND METHODS 

We arrange infinitely long filaments in the plane of a two-dimensional triangular lattice. 
The filaments are given an extensional spring constant a, and a filament bending modulus 
k. We introduce finite filament length L into the system by cutting bonds with probability 
1 — p, where < p < 1, with no spatial correlations between these cutting points. The 
cutting generates a disordered network with a broad distribution of filament lengths. When 
two filaments intersect, there exists a freely-rotating crosslink preventing the two filaments 
from sliding with respect to one another. Next, we introduce angular springs with strength 
K nc between filaments crossing at 60° angles with a probability p nc , where nc denotes non- 
collinear. These angular springs model the second type of crosslinker. See FigJT] for a 
schematic. 

We study the mechanical response of this disordered network under an externally applied 
strain in the linear response regime. For simplicity we set the rest length of the bonds to 
unity. Let be the unit vector along bonds and Uij = Ui — Uj the strain on the bond ij. 
For small deformation u, the deformation energy is 



where Pij is the probability that a bond is occupied, ^2uj) represents sum over all bonds and 
Y2(ijk) represents sum over pairs of bonds sharing a node. The first term in the deformation 
energy corresponds to the cost of extension or compression of the bonds, the second term to 
the penalty for the bending of filament segments made of pairs of adjacent collinear bonds, 
and the last term to the energy cost of change in the angles between crossing filaments that 
meet at 60° angle. Furthermore, for small deformations A6ijk = (iiji x rji — Ujk x rjk)-(fji x 
Tjk) = — ( M J»- r V+ M jfc- y V) _|_ Uik .r ik . It is straightforward to see that the angular spring ijk 
between ij and jk will contribute to an effective spring in parallel with ik, giving rise to an 
enhanced effective spring constant ji = a + §K nc . 




(ijk=n) 




(1) 



(ijk=ir/3) 
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A. Effective medium theory 



We study the effective medium mechanical response for such disordered networks follow- 
ing the mean field theory developed in \v\ for central force networks and [25] for filament 
bending networks. The aim of the theory is to construct an effective medium, or ordered 
network, that has the same mechanical response to a given deformation field as the depleted 
network under consideration. The effective elastic constants are determined by requiring 
that strain fluctuations produced in the original, ordered network by randomly cutting fila- 
ments and removing angular springs vanish when averaged over the entire network. 

Let us consider an ordered network with each bond having a spring constant fi m , a filament 
bending constant for adjacent collinear bond pairs K m , and an angular bending constant 
Knc,m between bonds making 60° angles. Under small applied strain, the filament stretching 
and filament bending modes are orthogonal, with stretching forces contributing only to 
deformations along filaments (un) and bending forces contributing only to deformations 
perpendicular to filaments (u±), and hence we can treat them separately. The angular 
forces due to the angular (non-collinear) springs, when present, contribute to stretching 
of filaments as discussed earlier, where we only consider three body interactions. For these 
springs to contribute to bending one needs to consider four-body interactions which is outside 
the scope of this paper and will be addressed in future work. 

We start with the deformed network and replace a pair of adjacent collinear bonds with 
bending rigidity K m by one with a rigidity k, and a bond spring with extensional elastic 
constant // m by a spring with an elastic constant /i and the facing 60° angular spring by 
K nc . This will lead to additional deformation of the above filament segments and the angle 
which we calculate as follows. The virtual force that needs to be applied to restore the 
nodes to their original positions before the replacement of the bonds will have a stretching, 
a bending and an angular contribution: F s , F&, and Fg. The virtual stretching force is given 
by F s = (fi m — a — 3K nc /2)u^ m , the virtual filament bending force is F^ = (K m — K)ux, m , while 
the virtual force to restore the angle is Fg = (K nc ,m — n nc )9 m , where u\\ >m , Ux,m an d 9 m are 
the corresponding deformations in the ordered network under the applied deformation field. 
By the superposition principle, the strain fluctuations introduced by replacing the above 
bending hinges and bonds in the strained network are the same as the extra deformations 
that result when we apply the above virtual forces on respective hinges and segments in the 
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unstrained network. The components of this "fluctuation" are, therefore, given by: 



d£n = 



IJ m /a* - fj, m + a + (3/2) 



K m /b* 



d6 = TTZ^ ( 2 ) 

The effective medium spring and bending constants, /i m , n m and « nc ,m, respectively, can 
be calculated by demanding that the disordered-averaged deformations (d£\\), (d£\\), and (dd) 
vanish, i.e. / — ^"""frf ,o ) = °> ( — £ m ~" z ) = 0, and / — \ = q. To 

\fJ,m/a*-fl m +a+3K, nc /2/ ' \ Km/» *-K m + K / ' \ Knc,m/c*—lt„- -->-"■-- I 



'"■n.c. rn T^nc 



perform the disorder averaging, since the stretching of filaments is defined in terms of spring 
elasticity of single bonds a, the disorder in filament stretching is given by P{a!) = p5(a'—a) + 
(1 — p)8(a'). Filament bending, however, is defined on pairs of adjacent collinear bonds with 
the normalized probability distribution P(k') = p 2 5(k' — k) + (1 — p 2 )5(k'). Similarly, for the 
angular springs, the normalized probability distribution is given by P(K,' nc ) = p nc p 
K nc ) + {\—p nc p 2 )5(K' nc )). This disorder averaging gives the effective medium elastic constants 
as a function of p and p nc as 

j ( „ -a-3n /2 \ + (J _ ( , m - 3^/7 



Hm/ a* \i m + a + 3/t arp /2 y \ji m / a* ji m + 3/t arp /2 

^rn P b ^-q^J ^"m,arp Pctrp P C 

K> 1 b* ^arp 1 ^* 

The constants a*, b* and c* for the network contribution to the effective spring con- 
stant Urn/ a* of bonds, to the filament bending rigidity K m /b*, and the bending rigid- 
ity n nc /c* of angular springs making 60° angles respectively, are given by a*,b*,c* = 
l^J2qTr [D Sjb>nc (q)D' 1 (q)~\ . The sum is over the first Brillouin zone and z is the coordi- 
nation number. The stretching, filament bending and non-collinear bending contributions, 
D s ^ nc {q) respectively, to the full dynamical matrix D(q) = D s (q) + D b (q) + D nc (q), are 
given by: 
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D s (q) = fi m ^2[l-e iqTii ] r tj r tj 

(ij) 

Db{q) = K m ^2 ~~ cos (9- r y)) 

(ij) 

-(1 - cas(2q.rij))] (I - r ijrij ) 
g 

D nc {q) = -K nCjm 2j [2(1 - cos(q.r^)) + 2(1 - cos(qr.r ifc )) 

-2(1 - cos(g.r jfe ))] ryr tt (4) 



with i" the unit tensor and the sums are over nearest neighbors 16) . Note that for small 
q, Db ~ q A and D s ~ q 2 have the expected wavenumber dependencies for bending and 
stretching. 

By definition, a* + b* + c* = 2d/z, where d = 2 is the dimensionality of the system. At the 
rigidity percolation threshold p = p rp , /i m , n m and « nC|TO vanish, giving a* = p+p 2 p n c—p 3 Pnc, 
b* = p 2 and c* = p 2 p n c- For semiflexible filament networks with only freely-rotating crosslinks 
i.e. filament stretching and bending interactions only, the rigidity percolation threshold is 
given by p rp = 0.457. For networks with angle-constraining crosslinks, at p nc = 1, we 
obtain rigidity percolation thresholds p rp = 0.405 for the case of flexible filament networks, 
and p rp = 0.347 for semiflexible filament networks. We also calculate how p rp changes on 
continuously increasing p nc from to 1. 

B. Numerical Simulations 

Simulations were carried out on a triangular lattice with half periodic boundary conditions 
along the shear direction for the energetic terms whose small deformation limit is given in 
Eq. ([1]). Networks were constructed by adding bonds between lattice sites with probability 
p. Next, a shear deformation was applied to the two fixed boundaries of magnitude ±7. 
The lattice was then relaxed by minimizing its energy using the conjugate gradient method 



301 ] allowing the deformation to propagate into the bulk of the lattice. Once the minimized 
energetic state was found within the tolerance specified, in this case the square root of 
the machine precision ~ 10~ 8 , the shear modulus was then measured using the relation, 
G = 2j5 7"ns , using small strains < 5%, with L denoting the system length and a ce u denoting 
the area of the unit cell for a triangular lattice which is equal to ?>\/2 in our units. System 



size L = 64 was studied, unless otherwise specified, and sufficient averaging was performed. 



III. RESULTS 



Mechanical integrity as measured by the shear modulus: On a triangular lat- 



tice, networks made sole 
around p rp j = 2/3 
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Hookean springs lose rigidity at a bond occupation probability 
321 ] . This result corresponds to the central force isostatic point 
at which the number of constraints is equal to the number of degrees of freedom on average. 
In contrast, networks made of semiflexible filaments become rigid at a smaller p due to ex- 
tra constraints placed on the system via filament bending. For semiflexible networks with 
freely-rotating crosslinks, our effective medium theory shows that the shear modulus, G, 
approaches zero at p rp = 0.457 as shown in Figj2] (a). This result is in good agreement with 



our simulation results yielding p rp = 0.442(6) and previous numerical results [26|. See FigJ2] 



(d). A different formulation of the EMT yields p rp ~ 0.56 [26|]. By introducing additional 
crosslinks that constrain angles between filaments at 60°, the rigidity percolation thresh- 
old is lowered. Our EMT yields p rp = 0.347 and our simulations yield p pr = 0.348(4) for 
Pnc = 1 (Figj2] (c) and (/)). The cooperative mechanical interplay between these crosslinks 
and their interaction with filaments allows the network to form a rigid stress-bearing struc- 
ture at remarkably low crosslinking densities, almost immediately after it attains geometric 



percolation, p c = 2sin(7r/18), which agrees with a calculation by Kantor and Webman [33 ]. 
For flexible filament networks, introducing angle-constraining crosslinkers also lowers the 
rigidity percolation threshold as compared to the isostatic point with the network attain- 



ing rigidity at p rp = 0.405 for our EMT and p rp = 0.408(4) in the simulations 



((FigO (6) 

and (e)). Incidentally, our result agrees very well with a previous simulation 34j. We also 
compute analytically and numerically how p rp changes with p nc . See FigJ^a), (6) and (c). 
Note that p rp is lowered continuously as the concentration of angle-constraining crosslinks 
is increased. 

Just above the rigidity percolation threshold, for a semiflexible network with freely- 
rotating crosslinks, we find a bending-dominated regime for sparse networks with the shear 
modulus eventually crossing over to a stretch dominated affine regime at higher filament 
densities. The purely stretch dominated regime is represented by the macroscopic shear 
modulus G staying almost constant with increasing p, while in the purely bend dominated 



regime the network is highly floppy and G is a sensitive function o 
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p, decreasing rapidly as 
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26J. For k < a, 



p is lowered. This behavior has been observed previously in 
both the effective medium theory and the simulations yield a bend-stretch coupled regime, 
which is characterized by an inflection in G as a function of p as observed most clearly for 
K = 10 -6 (with a = 1). 

We find a similar non-affine to affine crossover for the compositely crosslinked flexible 
filament networks and semflexible filament networks as p is increased. For the flexible 
filament networks, however, the bend-stretch coupling regime occurs for n nc <C a, i.e. K nc 
replaces k. For semiflexible filament networks, as long as n nc < k << a, the bend-stretch 
coupled regime is robust (for fixed p nc ). In contrast, for k « n nc « a, the angle- 
constraining crosslinker suppresses the bend-stretch coupled regime and enhances the shear 
modulus to that of an affinely deforming network (for fixed p nc )- The mechanics of the 
network has been altered with the introduction of the second type of crosslinker. 

Non-affinity parameter: To further investigate how the interaction of the crosslinkers 
affects the affine and non-affme mechanical regimes, we numerically study a measure for the 



degree of non-affinity in the mechanical response, T, defined in Ref. 26] as: 



1 - 

r = ZI 7 2 E(^- u -//) 2 - ( 5 ) 

i 

The non-affinity parameter can be interpreted as a measure of the proximity to criticality, 
diverging at a critical point as we approach infinite system size. We find that T develops 
a peak at the rigidity percolation threshold, which progressively moves to smaller values of 
p as the concentration of angular crosslinkers p nc is increased (FigJH (a)). A second peak 
develops near the isostatic point for n nc < k << a as seen in Figj4](6). As both the collinear 
and non-collinear bending stiffnesses tend to zero, the network mechanics approaches that of 
a central force network, and the second peak in T at the isostatic point becomes increasingly 
more pronounced. 

On the other hand, this second peak can be suppressed by increasing n nc j k (FigJH (6)), 
or by increasing the concentration p nc (FigJH (a)) even for very small values of n/a. This 
further corroborates that adding angle-constraining crosslinkers to non-affine networks can 
suppress non-affine fluctuations, provided they energetically dominate over filament bending. 
The reason for this suppression can be understood by considering the effect of adding a 
constraint which prohibits the free rotation of crossing filaments. As the concentration of 



these non-collinear crosslinks p nc is increased (at fixed avg. filament length) microscopic 
deformations will become correlated. The lengthscale associated with this correlation will 
increase on increasing either p or p nc , and will eventually reach a lengthscale comparable to 
system size even at p ~ p T p,i at large enough concentration and/or stiffness of the angular 
springs. As a result the mechanical response of the network will approach that of an affinely 
deforming network. Upon decreasing the value of K nc /a relative to n/a we again recover 
the second peak because energetically the system can afford to bend collectively near the 
isostatic point. 

Scaling near the isostatic point: Finally, using scaling analysis we quantify the simi- 
larity in mechanics between freely-rotating crosslinked semiflexible networks and compositely 
crosslinked flexible networks. To do this, we examine the scaling of the shear modulus G 



near the isostatic point with Ap = p — p rp j <C 1. For k/ol <C Ap (or K nc /a <C Ap 



the 



shear modulus scales as G = a\Ap\fg±(£ |Ap| - *) (or G = a|Ap|^±(^|Ap|-*)) |26|, |35|. 
For both (a) k = 0,K nc > and (b) k > 0, n nc = 0, the EMT predicts / = 1 and <fi = 2 as 
shown in Figj5](a) and (b), indicating that both types of networks demonstrate redundant, or 
generic, mechanics. To compare the EMT results with the simulations, we use the position 
in the second peak in F to determine the central force percolation threshold, p rp ,i, and then 
vary / and to obtain the best scaling collapse. For CSjSG (&) i Prp,I 

= 0.666(3), / = 1.1(1) 

and <p = 2.8(1). For case (b), p rp>1 = 0.659(5), / = 1.1(1) and <p = 2.9(1). Both sets of expo- 
nents are reasonably consistent with those found in Ref. (26 1 for a semiflexible network with 
freely-rotating crosslinks only. Preliminary simulations for compositely crosslinked semiflex- 
ible networks indicate that the shear modulus scales as G = a\ Ap\fQ±(^\ Ap\~^, -^|Ap|~ 7 ) 
also with a similar / and a similar <fi with = 7. 



IV. DISCUSSION 



In the limit of small strain, we conclude that the presence of multiple crosslinkers in living 
cells can be simultaneously cooperative and redundant in response to mechanical cues, with 
important implications for cell mechanics. Redundant functionality helps the cytokeleton be 
robust to a wide range of mechanical cues. On the other hand, different crosslinkers can also 
act cooperatively allowing the system to vary the critical filament concentration above which 
the cytoskeleton can transmit mechanical forces. This may enable the cytoskeleton to easily 
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remodel in response to mechanical cues via the binding/unbinding of crosslinkers (tuning 
concentration) or their folding/unfolding (tuning stiffness and type of crosslinker) . Since the 
cytoskeleton consists of a finite amount of material, the ability to alter mechanics without 
introducing major morphological changes or motifs may play important role in processes 
such as cell motility and shape change. 

Cooperativity: In our study of two types of crosslinkers, crosslinkers that allow free 
rotations of filaments and crosslinkers that do not, we find two types of cooperative effects 
in the mechanics of such compositely crosslinked networks. The first cooperative effect de- 
pends on the relative concentration of the two types of crosslinkers and second depends on 
the relative stiffness of the angle-constraining crosslinkers to the bending stiffness of the 
individual filaments. The first cooperative effect can be most strikingly observed beginning 
with an actin/a-actinin network and increasing the concentration of FLNa, with a-actinin 
representingthe freely-rotating crosslinker [6] and FLNa representing the angle-constraining 
crosslinker [7|. By tuning the concentration of FLNa, the cell can modulate the minimum 
concentration of actin filaments necessary to attain mechanical rigidity, which can be es- 
sentially as low as the filament concentration required to form a geometrically percolating 
structure. This is in good agreement with the experimental observation that FLNa creates 
an F-actin network at filament concentrations lower than any other known crosslinker 
Increasing the FLNa concentration also suppresses the non-affine fluctuations near the rigid- 
ity percolation threshold by increasing the shear modulus of the network and giving rise to 
a more affine mechanical response while keeping the filament concentration fixed. Moreover, 
the cooperativity of a-actinin and FLNa working to ehance the mechanical stiffness of actin 
networks has recently been observed in experiments 27] . The addition of angle-constraining 
crosslinkers to flexible filament networks also decreases the concentration threshold required 
for mechanical rigidity, though the lower bound on the threshold is not as close as to geo- 
metric percolation as it is for semiflexible filaments. The lowering of the rigidity percolation 
threshold is independent of the energy scale of the crosslinker. It depends purely on the 
number of degrees of freedom the crosslinker can freeze out between two filaments, i.e. the 
structure of the crosslinker. 

The second cooperative interplay between the two crosslinkers depends on the energy scale 
of the angle-constraining crosslinker to the filament bending energy. For kCch, the freely- 
rotating semiflexible filament system exhibits large non-affine fluctuations near the isostatic 
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point. Upon addition of the angle-constraining crosslinkers, for K nc > k, the non-affine 
fluctuations near this point become suppressed and the mechanics of the angle-constraining 
crosslinker dominates the system. Once again, with a small change in concentration of the 
second crosslinker, the mechanical response of the network is changed dramatically. 

Redundancy: We observe two redundant effects in these compositely crosslinked net- 
works, the first of which depends on energy scales. For n nc <C k with nCa, the non-affine 
fluctuations near the isostatic point in the freely-rotating crosslinker semiflexible filament 
network remain large even with the addition of the angle-constraining crosslinker. In other 
words, the angle-constraining crosslinkers are redundant near the isostatic point. Their pur- 
pose is to decrease the amount of material needed for mechanical rigidity as opposed to alter 
mechanical properties at higher filament concentrations. 

Redundancy is also evident in the mechanics of these networks sharing some important, 
generic properties. All three networks studied here (free-rotating crosslinked semibflex- 
ible networks and compositely crosslinked semiflexible and flexible networks) have three 
distinct mechanical regimes: a regime dominated by the stretching elasticity of filaments, a 
regime dominated by the bending elasticity of filaments and/or stiffness of angle-constraining 
crosslinkers, and an intermediate regime which depends on the interplay between these in- 
teractions. The extent of these regimes can be controlled by tuning the relative strength of 
the above mechanical interactions. In particular, the ratio of bending rigidity to extensional 

n 

modulus of an individual actin filament is ~ 10~ 3 [21]. Since the bend-stretch coupled 
regime has not been observed in prior experiments on in-vitro actin networks crosslinked 
with FLNa only, we conjecture that the energy cost of deformation of angles between fila- 
ments crosslinked with FLNa is larger than the bending energy of filaments. The qualitative 
redundancy becomes quantitative, for example, near the isostatic point where we obtain 
the same scaling exponents for G as a function of p — p rp j and k(ot K nc ) for the free- 
rotating crosslinked semiflexible network and the compositely crosslinked flexible network. 
Preliminary data suggests the same scaling extends to compositely crosslinked semiflexible 
networks. This result is an indication of the robustness of these networks and should not be 
considered as a weakness. Whether or not this robustness extends to systems experiencing 
higher strains such that nonlinearities emerge is not yet known. 

Lamellipodia mechanics: The interplay between cooperative and redundant mechani- 
cal properties may be particularly important for the mechanics of branched F-actin networks 
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in lamellipodia. Within lamellipodia, there exist some filament branches occuring at an an- 
gle of around 70° with respect to the plus end of the mother filament (referred to as Y — 

n 

junctions). These branches are due to the ABP Arp2/3 [101 ] . During lamellipodia formation, 
these branches are presumed to be the dominant channel for filament nucleation. The me- 
chanics of Arp2/3 can be modeled as an angular spring between the mother and daughter 
filament with an angular spring constant of approximately 10~ 19 J rad~ 2 lOj. In other words, 
Arp2/3 is an angle-constraining crosslinker for Y — junctions (as opposed to X— junctions), 
and thereby plays an important role in lamellipodia mechanics as demonstrated in this work. 
The mechanical role of Arp2/3 in lamellipodia has not been inve stig ated previously and may 
help to discriminate between the dendritic nucleation model [ll[|l2| and a new model 13] by 
predicting the force transmitted in lamellipodia as a function of the Arp2/3 concentration. 

In addition to Arp2/3, FLNa localizes at A— junctions in the lamellipodia and is thought 
to stabilize the dendritic network ]36j. Both angle-constraining crosslinkers lower the fila- 
ment concentration threshold required for mechanical rigidity in the system. Depending on 
the energy scale of FLNa as compared to the energy scale of Arp2/3, addition of the FLNa 
may or may not modulate, for example, the bend-stretch coupling regime at intermediate 
filament concentrations. Again, at times mechanical redundancy is needed and at times not. 
With three crosslinkers, the system can maximize the redundancy and the cooperativity. Of 
course, lamellipodia are dynamic in nature and are anisotropic since the Arp2/3 is activated 
from the leading edge of a cell. Both attributes will modulate the mechanical response. 

Outlook: We have demonstrated both cooperativity and redundancy in the mechanics 
of compositely crosslinked filamentous networks. We have done so while maintaining the 
structure of an isotropic, unbundled filament network. Of course, crosslinkers can alter the 
morphology of the network via bundling, for example. In other words, different crosslinkers 
serve specific functions. This specificity results in a change in microstructure. This will pre- 
sumably affect the mechanics such that the cooperative and redundant interactions between 
multiple crosslinkers may differ from the above analysis. For example, the crosslinker that 
dominates in terms of creating the morphology will presumably dominate the mechanics. 
Schmoller and collaborators [5| suggest that crosslinker with the higher concentration deter- 
mines the structure and, therefore, the mechanics. Instead of redundancy or cooperativity, 
the specificity leads to the simple additivity of two types of crosslinkers in that different 
crosslinkers act independently of one another. In this study, however, we find both coop- 
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erativity and redundancy in the network mechanics even in the absence of such structural 
changes [37j, which, is arguably less intuitive and, therefore, more remarkable. Finally, while 
our focus here has been on the actin cytoskeleton as an example of a filamentous network, 
our results can be extended to collagen networks as well 
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FIG. 1: Deformed configuration a compositely crosslinked semiflexible network with 2.7 percent 
strain, with bond occupation probability p = 0.64, and angle-constraining crosslinker occupation 
probability p nc = 0.15 The purple lines denote semiflexible filaments, the red arcs denote angle- 
constraining crosslinks, the black circles represent nodes where all crossing filaments are free to 
rotate, while the grey circles denote nodes where some of the crossing filaments are free to rotate. 
The filament bending stiffness relative to stretching stiffness k/ol = 10~ 6 and the stiffness of angular 
crosslinks relative to stretching stiffness K nc /a = 10 -6 . 
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FIG. 2: The shear modulus as a function of p for semiflexible networks with freely-rotating 
crosslinks ((a) and (d)), flexible networks with freely-rotating and angle-constraning crosslinks 
((b) and (e)), and semiflexible networks with both crosslinkers ((c) and (f)).The top panels show 
results from the effective medium theory and bottom panels show results from the simulations. 
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FIG. 3: The presence of angular constraints allows these networks to have a finite rigidity even for 
small concentration of filaments. Figure (a) shows how the rigidity percolation threshold can be 
continuously lowered by increasing the concentration of angular springs for flexible (solid, blue) and 
stiff (dashed, red) networks. The lines correspond to the effective medium theory and the symbols 
to the numerical simulation, where the system size L = 32 for the semiflexible filaments and L = 64 
for the flexible filaments. Figures (b) and (c) show the shear modulus (in logarithmic scale described 
by the colorbar) as a function of p and p nc for flexible networks (b) and semiflexible networks (c). 
The parameter values studied are (b) K nc /a = 10~ 4 and (c) k/ol = 10~ 4 , K nc /a = 10~ 2 . The black 
dashed lines in (b) and (c) correspond to the effective medium theory prediction of the rigidity 
percolation threshold. 
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FIG. 4: The non-affinity parameter T as a function of p for semiflexible networks with both types of 
crosslinkers. In (a) we show the effect of changing the concentration p nc of the angle-constraining 
crosslinkers for n/a = 10 -4 , K, n c/d = 10 -2 and L = 32, while in (b) we show the effect of changing 
their stiffness K nc for L = 64. 
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FIG. 5: Close to isostaticity, the shear modulus G scales with Ap = p — p rp j and n (fv nc ) as 
G\Ap\~f = K,\Ap\~^. The effective medium theory predicts mean field exponents / = 1 and <f> = 2 
for both semiflexible networks with freely-rotating crosslinkers (a) and compositely crosslinked 
flexible networks (b), while simulations predict / = 1.1(1) and (p = 2.9(1) for semiflexible networks 
with freely-rotating crosslinkers (c) and / = 1.1(1) and <ft = 2.8(1) for compositely crosslinked 
flexible networks (d). 
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